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Abstract 



We consider the classical equations of motion of SU(2) gauge theory, with- 
out a Higgs field, in Minkowski space. We work in the spherical ansatz and 
develop a perturbative expansion in the coupling constant g for solutions 
which in the far past look like freely propagating spherical shells. The topo- 
logical charge Q of these solutions is typically non-integer. We then show that 
Q can be expressed as a power series expansion in g which can be nonzero 
at finite order. We give an explicit analytic calculation of the order g 5 con- 
tribution to Q for specific initial pulses. We discuss the relation between 
our findings and anomalous fermion number violation, and speculate on the 
physical implications of our results. 
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I. INTRODUCTION 



Fermion number is not conserved in the standard model U, yet reliable techniques for 
calculating fermion number violating rates in high energy collisions have not been fully 
developed. At very low energy where the processes are best described as tunnelling events, 
the calculational methods Q make use of solutions to the Euclidean field equations - 
instantons [0. At energies comparable to but below the sphaleron barrier ||, Euclidean 
methods [|J and other methods in which part of the calculation is done in Euclidean space 
H have also been applied. However, at energies well above the tunnelling barrier, it may be 
more appropriate to work directly in Minkowski space ||[7] . To this end, we are investigating 
solutions to the Minkowski space classical equations of motion. In this paper, we study SU (2) 
gauge theory without the Higgs field. 

Working in the (spatial) spherical ansatz, we explore solutions which have the property 
that in the far past they describe freely propagating incoming shells of energy. We do this by 
developing a perturbative expansion in the coupling g which can be used to systematically 
solve the equations of motion once arbitrary initial profiles have been specified. We discuss 
the topological charge of these solutions, which we also develop in a power series expansion 
in g. We show that the topological charge is nonzero at order g 5 . We then investigate the 
associated anomalous fermion production in the presence of these solutions. 

Recently , it was shown that in the spherical ansatz the classical equations of motion of 
SU (2) gauge theory can be reduced to two coupled nonlinear partial differential equations of 
motion for two gauge invariant functions of r = |x| and t. In Ref. J7[, new solutions to these 
equations were also presented which have much in common with the previously discovered 
solutions U of Liischer and Schechter (LS). As we will show here these explicit solutions are 
examples of a wide class of finite energy solutions all of which have certain general features 
in common. At early times they depict a thin spherical shell of energy imploding towards 
the origin at near the speed of light. At intermediate times the region around the origin 
is energetically excited and at late times the shell is expanding outward, asymptotically 
approaching the speed of light. 
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We are particularly interested in solutions which look free in the far past. This is because 
our ultimate goal is to connect the classical solutions to the quantum mechanical description 
of the scattering of quanta which in the far past are well separated and non-interacting. For 
a classical solution to be said to be free in the far past we require that there exists a gauge 
in which the commutator term in (see 1.3) can be neglected in the far past. This means 
that the theory linearizes. We also want A v in this gauge to be independent of g in the far 
past. This means that the initial value data can be specified without reference to g. The 
LS solutions and the solutions of Ref. 0] do not satisfy this requirement. However, in this 
paper we will show, in the spherical ansatz, how to specify initial data which do meet this 
criterion. 

Returning to the description of the shells, if we focus on the region near the origin we see 
that before the front reaches it, it is devoid of energy, i.e. it is a vacuum configuration. The 
energy shell excites this region and then leaves so that at late times the region is once again 
a vacuum configuration. However, the late time vacuum structure of this region need not 
coincide with its early time structure. As we now explain, the region can be associated with 
a local topological charge which is typically nonvanishing, non- integer, and ^-dependent. 

We define the local topological charge inside a sphere of radius R as 

Q(R) = -f^ r dt J d 3 xe^ at3 Tr(F^F a(3 ) , (1.1) 

3271^ J-oo J R 

where the spatial integration is over the interior of the sphere. For the LS solutions, Q(R) 
approaches a fixed value as R goes to infinity whereas for the solutions of Ref. [0], Q(R) 
grows without bound. Let us restrict our attention to all those solutions for which Q(R) 
has a limit as R goes to infinity. In this R goes to infinity, Q(R) approaches the 

topological charge Q defined by 

Q = 4-2l d ' x ^Tr(F M ,F a/3 ) , (1.2) 

where the integration is over all of space-time. 

The topological charge Q can take any value, not just integer values [0. At first sight 
this may appear to contradict the observation that (|1.2p is a topological invariant of the field 
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configuration. Let us review the usual argument || that leads to integer topological charge. 
Imagine that the region of space-time where the energy density is nonzero is bounded. 
Surround this region with a three dimensional surface which is topologically S 3 . On this 
surface = so the gauge field is pure gauge, i.e. = (i/g) Ud^W where U takes values 
in 577(2). This defines a map from S 3 into 577(2) which is characterized by an integer since 



773(577(2)) = Z. The integrand in (1.2) is a total divergence so Q can be written as an integral 
over the surrounding three dimensional surface. This integral computes the winding number 
of the map from S 3 into 577(2) and therefore Q is an integer. 

This argument assumes that the space-time region where the energy density is nonzero 
is bounded. For solutions to the equations of motion, energy is conserved and the energy 
computed on any equal time surface is non-vanishing. It is therefore impossible to surround 
the space-time region of nonzero energy density with a three sphere. We have no reason to 
expect and we do not find integer values of Q. 

If the topological argument applied and Q were an integer then certainly Q would be 
independent of the gauge coupling g. We now examine another commonly used argument 
which, when applicable, says that Q is independent of g. The field strength F^ v is defined 
as 

F„ v = d^A u - d u A^ - ig[A^ A u ] , (1.3) 
and the equation of motion is 

D»F, U = d»F, u - ig[A», F^] = . (1.4) 

We can redefine the fields so that the equation of motion is independent of g. Let 

A, = -A', , (1.5) 
9 

F^v = -F\ v = - (d,A' v - d v A\ - i[A^, A',]) . (1.6) 
Now the equation of motion reads 



& t F J IJU/ -i[A" l ,F' lw ] = 0, (1.7) 

which contains no g. The topological charge (|1.2|) is seen to be independent of g: 

Q[A'»] = 3-^5 / d'xe^Tr (F'^F*^) . (1.8) 

In this sense solutions to the equation of motion yield a topological charge independent of 
the coupling. 

However, we are interested in solving (|1.4j) in Minkowski space with fixed initial value 
data which does not depend on the coupling g. That is, we wish to solve ( |1.4| ) specifying 
the field and its time derivative A^ at some fixed time T. We specify 

A M (x,T) = Aj*(x), (1.9) 
i,(x,T) = V; nit (x) , (1.10) 

where y4™ [t (x) and V^ mt (x) are arbitrary functions consistent with the constraints. In terms 
of rescaled variables 

AV(x,T) = ^jf(x) , (1.11) 

i' M (x,T) = ^; nit (x). (1.12) 

If we fix A^^x) and V^ mt (x) then A ^ actually depends on g. Accordingly Q depends on g. 

As an example, if g is equal to zero, we are solving ( |1 . 7|) with the condition that A ^ and 
A ^ are zero at T. The solution is A' M (x, t) = which yields zero topological charge in fll.8| ). 
Another way to see this is to work with g = directly. The equation of motion (|1.3|) is 



d»(d ll A v -d v A li ) = 0, (1.13) 

which we can solve with fixed initial value data, yielding non-trivial solutions. In this case 
Q given by ( |1.2| ) is seen to be zero since g is zero. Thus if we solve the classical equations 
of motion with fixed, (/-independent, initial value data we expect the topological charge to 
vanish as g goes to zero. This suggests that we can develop a power series expansion in g 
for Q. 



To conclude this introduction, we describe the organization of the rest of the paper. In 
section II we review the spherical ansatz and write the equations of motion in terms of two 
gauge invariant variables. In section III we show that in the spherical ansatz, the theory 
linearizes in the far past and the far future. We express the requirement that A^ 1 be g- 
independent in the far past in terms of the gauge invariant variables. In section IV we show 
how to solve the field equations as a perturbative expansion in g for arbitrary initial pulses, 
and perform the calculation explicitly to zeroth and first order in g. In section V we return 
to our discussion of the topological charge of the solutions. We show how to expand Q as a 
power series in g. We present examples of initial data for which the topological charges of 
the corresponding solutions are nonzero at order g 5 . In appendix A, without restricting to 
the spherical ansatz, we show that in source-free electromagnetism the topological charge of 
localized pulses vanishes. In section VI we discuss anomalous fermion number violation in 
the presence of classical gauge field sources. We show that the number of fermions produced, 
in a quantum average sense, is Q. In section VII, we remark upon the implications of our 
results for anomalous fermion number violation in high energy scattering in the Standard 
Model, and discuss open questions. 

II. THE SPHERICAL ANSATZ 



Here we review the spherical ansatz JT0| for SU (2) gauge theory without a Higgs field and 
show how the equations can be reduced to two equations for two gauge invariant functions. 
The spherical ansatz is given by expressing the gauge field in terms of four functions 
ao , a\ , a and 7 of r and t: 

A (-x,t) = ~a (r,t)<r ■ x 

1 otiv^ t\ 7(7*, t\ 
Ai(-x,t) = ~[ai(r,t)(T ■ H — - a ■ xx;) H — e ijk Xja k ] , (2.1) 

where x is the unit three- vector in the radial direction. Note that we do not introduce 1/g 
factors in (12.11) as was done in Refs. [BHIl and [ffll. The action 
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S = ~Jd 4 xTr(F^Fn 



(2.2) 



in the spherical ansatz takes the form of an Abelian Higgs model in curved space |T0 | . 

2" 



S = 4vr / dt j dr 



where 
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1 



2r 2 \ g 2 



(2.3) 



(2.4) 



is the field strength, 



is a complex scalar and 



X = a + i(j- -) , 



(2.5) 



(2.6) 



is the covariant derivative. The indices are raised and lowered with the 1 + 1 dimensional 
metric ds 2 = —dt 2 + dr 2 . (Note that in this paper we raise and lower indices with a flat 
space metric; we do not absorb the explicit factors of r into a curved space metric as was 
done in Ref. 0.) 

The equations of motion for the reduced theory are 



d»{r 2 f, v )=ig [(D vX Tx-X*DuX\ 



(2.7a) 



and 



D 2 + - 2 (g 2 \x\ 2 - i) 



x = o 



It is convenient to write the complex field x i n polar form 



(2.7b) 



X(r,t) = -ip(r,t)exp[ig<p(r,t)] , 



(2.8) 



which casts the equations of motion in terms of p , cp and Qj^ clS 



& t (r z f lu ,) + 2gY(d v <p-a„) = 



^2.9a) 



d^p - g 2 p(d^ - a")(0 M <p - a„) - -^V - l)p = 



(2.9b) 



and 



5" 



p 2 {d^ - a,) 



. 



:2.9c) 



Note that (|2.9cj ) follows from (|2.9a| ) so there are three, not four, independent equations of 
motion which is expected because of the residual U(l) gauge invariance. 

In 1 + 1 dimensions f^ u must be proportional to e Mi , so we can define the gauge invariant 
field if) through 



r 2 f^u = -2e Mi/ V 



(2.10) 



(Here e i = +1.) Equation ( |2.9a|) now becomes 



d a if, = -g* e au p 2 (d u <p - a 



(2.11) 



which implies 



p 2 



Substituting ( PH) into (|2~9bl) gives 



1 



^p + — ^ ^ - -o?y - i)p = o 



(2.12a) 



(2.12b) 



Equations ( |2.12a| ) and ( |2.12b|) are the equations of motion for the gauge invariant vari- 
ables p and if). If a gauge is chosen then the gauge-variant variables <p and can be 
determined using (|2.11|) . For example, in the <p = gauge, given p and if>, equation ([2.11 



directly determines a v . Equations (|2.12a| ) and ( |2.12b| ) for if) and p can be solved after spec- 
ifying initial value data, that is the values of if> , if> , p and p at some fixed time. Any initial 
value data expressed in terms of if) , if) , p and p will be consistent with the Gauss's law 



S 



constraint. This is because Gauss's law is the v = component of equation ( |2.9a| ) which is 
equivalent to the a = 1 component of ( [2. lip . 

Using the equations of motion, the energy associated with the action ( |2.2| ) can be written 
in terms of p and ip as 



E = 8tt / dr 
Jo 



(dtp) 2 + ^{drpy + jjWW + tt^Ow + Sr + 



V 2 , QrV-i) 5 



2^ p 



Ag 2 r 2 



(2.13) 



From this expression we see that the field configuration with zero energy has ?/> = and 
P = 1/9- 

It is convenient to introduce a shifted field 5 defined by 



2gS=\-gV . 



(2.14) 



Equations (2.12) now become 



dM--ilj = -2gdu 



.1-20*, 



and 



^-^ 5 = ^T32oJ-^T32^-^ 5 ' 



and the energy ( |2.13| ) is 



E = 8vr / dr 
/o 



(2.15a) 



(2.15b) 



roc I , >, qlj 2 R 2 ' 



(2.16) 



.2(1-2^) 

We will show how to solve ( 2.15|) in a perturbative expansion in g after we show that 
these equations generically describe imploding shells in the far past which turn into outgoing 
shells in the far future. 
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III. LINEARIZATION IN THE FAR PAST AND FAR FUTURE 



In the introduction we pointed out that in the spherical ansatz a typical solution to the 
equations of motion describes an imploding shell, moving close to the speed of light, at very 



early times. To see how this comes about consider equations ( 2.15 ) and imagine that at 
some early time T < 0, 5 and ip are both pulses of width A centered at r near \T\ with 
A -C |T|. By a pulse we mean a function which is very close to zero except in a region whose 
size is A. Accordingly, space and time derivatives of ip and 5 are pulses with amplitude of 
order ip/A and 5/ A. For r ~ \T\ > A we can now neglect the 1/r 2 terms in ( |2.15| ). We 
then see that if ip(r,t) and 5(r, t) depend only on r + t, that is, i[)(r,t) = ip p (r + t) and 
5(r, t) = 5 p {r + i) then equations (|2.15|) are satisfied. If ip p {u) and 5 p {u) are close to zero 
except for u in a region of size A around u = 0, then the solution describes an incoming 
shell of width A moving undistorted along r = —t. This description remains valid for all 
t < -A. 

Although equations (|2.15|) with the 1/r 2 terms neglected have solutions ip(r, t) = ip p (r+t) 
and S(r,t) = 5 p (r + t), these equations are not linear. We are for example not assuming 
that 5 and ip are small. There is however a sense in which we can say that the theory 
linearizes. Consider the field strength in the 3 + 1 dimensional theory given by ( |1.3p . 
We now argue that there is a choice of gauge such that the commutator term, [A^-Aj,], 
can be neglected relative to the linear term, d il A v — d u A fl , for all t much less than —A. 
Clearly this statement requires a gauge choice since the size of and its derivatives are 
gauge dependent. Nonetheless the statement that there exists a gauge with this property is 
a gauge invariant characterization of solutions in the spherical ansatz. 

Work in the A = gauge which in the reduced theory is equivalent to the a = gauge. 
Using (|2.11| ) and (|2.14|) we then have 



d r il>(r,t) = (2g6{r,t)-l)dt<p{r,t) , (3.1) 



or 
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where <po(r) is an arbitrary function of r. Now for t <C —A, using ip p (r + t) and S p (r + t) 
for ip and 5, we have 

W*) = -/ 9 r 7 N +yo(0- 3 - 3 

J-oo 1 — 2gdp(u) 
Equations (|2.11 ) and ( 2.14|) also give, for all time, 



which for t —A, upon using (|3.3| ), becomes 

Oi(r,*) = <9 r v9 (r) . (3.5) 

In the ao = gauge we are still free to make a time independent gauge transformation which 
we do to set 

ai {r,t) = (3.6) 

for t —A. What we have shown is that in the ao = gauge if we set ai = at one early 
time then we have ai = at all early times. 

If we now turn to (|2 . 1|) we see that in the gauge we have chosen, at early times, the 
dominant pieces in Ai go like a/r and j/r. Since the gauge has been specified, a and 7 are 
determined by 6 and if) so, for example, derivatives of a and 7 are of order a/ A and 7/ A. 
This means that d^A^ is of order ct/(rA) and j/(rA) whereas A^Aj is of order a 2 /r 2 , ofj/r 2 
and 7 2 /r 2 . For A <C r the commutator term in can be neglected relative to the linear 
term. (In fact F^ v goes like l/(rA) so the total energy in the shell scales as 1/A.) We also 
see that the energy is, at early times, quadratic in the field A4 which is a signature of a 
linear theory. 

We are particularly interested in solutions which in the far past describe free particles. 
Thus we want the gauge field A^, in the gauge in which the solution linearizes in the far 
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past, to be independent of g in the far past. In this case is also independent of g since 
Ffu, ~ d IM A u — d u A fl . Then, for example, the energy which depends on is independent of 

9- 

In the gauge in which the solutions linearize, a® = and a\ = for t <C —A. For in 
this regime to be independent of g, we require that a and 7 be independent of g (see fl2.1|) ). 
Now from fl2j), Q and ( pll , we have 

5 = 7-^(a 2 + 7 2 ) . (3.7) 



Using O, O, ( gig) and Q, we also have 

= -d t a + g{jd t a - ada) . (3.8) 

This means that if in the far past a and 7 are independent of g, then the gauge invariant 

variables if> and 5 have pieces which are zeroth and first order in g. If a(r, t) and j(r, t) for 
t — A are described by pulses a(r, t) = a p {r + t) and 7(7-, t) = 7 P (r + t) then 

<5 P H = 7p(w) - 2 9 ( a p( u ) + % 2 ( M )) ( 3 - 9a ) 

and 

d u ip P {u) = -d u a p (u) + g{-f p (u)d u a p (u) - a p (u)d u >y p (u)^ . (3.9b) 



Note that if we assume this form for 5 p (u) and ip p {u), then the energy evaluated using ( p,16| ), 
neglecting the 1/r 2 terms, is indeed independent of g. 

We have shown that an incoming pulse moving undistorted at the speed of light, if)(r, t) = 
if) p (r + t) and <5(r, t) = 5 p {r + t), solves the equation of motion ( [2.15| ), as long as the pulse 
is located at a radius much larger than the width of the pulse. Similarly an outgoing pulse, 
if)(r,t) = ipf(t — r) and 5(r, t) = 5f(t — r) solves the equations of motion under the same 
conditions. A general scattering solution consists of an incoming pulse in the far past which 
in the far future is an outgoing pulse. The future pulse profiles, ipf and 5f will not, in 
general, coincide with the past pulse profiles ip p and 5 P . In fact a full understanding of the 
solutions to ( |2.15[) would allow us to determine if>f and Sf given arbitrary ip p and S p . 
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IV. THE PERTURB ATIVE EXPANSION 



In this section we solve equations (|2.15| ) perturbatively in g in an in-field formalism 



where we specify the form of the solution in the far past. If we set g equal to zero, equations 
fl2.15a| ) and (|2.15b| ) have the same form: 

-d 2 + d 2 ~^jij = 0, (4.1a) 

and 

2 



~ d t + 9 r ~ ^ ) 6 ° = • ( 41b ) 



It is convenient to switch to light cone coordinates defined as 

u = t + r 

v = t-r, (4.2) 

in which ds 2 = —dt 2 + dr 2 = —du dv and <9 M <9 M = —d 2 + d 2 = —4d u d v . Finite energy requires 
the solutions to vanish at r = (u — v)/2 — 0. The general finite energy solution to the g = 
homogeneous equations ( |i.l|) is 



Mu, v) = G'(u) + G'(v) — \G(u) - G(v)} (4.3a) 

u — V 

and 

5 (u, v) = H\u) + H\v) — [H(u) - H(v)) , (4.3b) 

u — V 

where G(x) and H(x) are two arbitrary functions which approach constants as x goes to 
±oo. If G'(x) and H'(x) are both nonzero only in a region of size A containing x = 0, 
then the solution represents an energy shell of width A moving along r = — t at early times 
(t -C —A) and along r = t at late times (t > A). 

We now expand the fields if) and 5 in the coupling g: 

^ = ^0 + ^1+^2 + ... , (4.4a) 
5 = 5 + g5 1 + g 2 5 2 + ... , (4.4b) 
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and substitute into the equations of motion Q2.15Q . Expanding the equations in g gives 

2 



d u d v + 
d u d v + 

d u d v + 
d u d v + 





-v) 2 _ 




2 


(u 


-v)\ 




2 


(u 


-v) 2 _ 




2 



^ = If 

Si = if 



£ TO 



(u — v) 2 

where 1% and l b n both depend only on ip , . . . ip n -i and Sq, . . . 8 n -\. For example 



(4.5) 



lf = ^(Sod^ ) 

if = i«Vo - \0^0 d^ + 

1\> _ 1 
~ 2' 



II = MSi d^ + 5 d^x + 25% d^o) 



(4.6a) 
(4.6b) 
(4.6c) 
(4.6d) 



2-"-"-"- 2 ~ u ru 2 ru ' 2~ / *" 1 ~ " u ' (m-w) 2 

Given the sources 1$ and equations ( |4.5| ) do not uniquely specify ^ n and 5 n because 
we can add homogeneous solutions of the form ( |4.3a| ) and ( |4.3b| ) to ip n and 5 n and still satisfy 
the equations. We will first solve ( |4.5| ) with the requirement that the solutions vanish in 
the far past, and will later use the freedom to add homogeneous solutions to find solutions 
which describe ^-independent a and 7 pulses in the far past. Consider the equation 

2 



d u d v 



(u — v) 



$(«, v) = I(u, v) 



(4.7) 



where the source I — > as t> — > — 00, which we solve requiring that the solution $ — > as 
v — > —00, and also that $ vanishes at r = (u — v)/2 = 0. The retarded Green's function 
solution is 



ru rv 

§(u,v)= du' dv' F{u,v]u',v')I{u',v') 

J V J — CO 



(4.8) 
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where 

(m — f — tr) 

(Note that F(u, v; u', v') satisfies Q) with / = 0.) 

We can use ( [4.8| ) to get explicit formulas for ip% and Si given the sources if and /f of 
( |4.6a| ) and ft4.6b|) along with ip and So expressed in terms of the arbitrary functions G(x) 
and H(x) of ( t4.3a|) and Q4.3b|) . The result is 



Mu,v) = -G\u)H\v) - G\v)H'{u) - . 6 .J G(u) - G(v))[H(u) - H(v)) 

[U — v) z 

+^—[G\u) + G'(v)][H(u) - H{v)\ 
u — V 

+-^[H'{u) + H'{v)][G{u) - G(v)) (4.10a) 



and 



Sfav) = G'{u)G'{v) + G .^ 2 - -*—[G\u) + G'(v)][G(u) - G(v)} 

[u — v y u — V 

+ JL_ [ U dx G'(x) 2 - H'(u)H'(v) - , 5 ,J H(u) - H(v)} 2 

U — V Jv {U — V) z 

+^—[H'(u)+H'(v)}[H(u)-H(v)} + ^— f U dxH'(x) 2 . (4.10b) 

U — V U — V Jv 

Using the fact that G'(x) and H'(x) are only nonvanishing near x = we see that ipi and 
Si do in fact go to zero as v — > — oo. We are still free to add homogeneous pieces to ipi and 
Si if we so choose. Indeed we will need to use this freedom since equations ( |3.9| ) and the 
requirement that a and 7 are (^-independent in the far past imply that ipo, So, and Si 
are in general nonzero as v — > —00. Given tpi and 8\ we can evaluate if and 1 "| in ( |4.6c|) 



and Q4.6d| ) and then obtain ip 2 and S2 using (|4.8|) . Continuing along these lines we can, in 
principle, get ip and S to any order in g. Of course in practice the calculations get more 
complicated with each order of perturbation theory. 
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V. TOPOLOGICAL CHARGE 



In the spherical ansatz, using the equations of motion, the local topological charge ( |1 . 1|) 
can be written [0] as 

Q(R) = £- r dt f R dr (-d 2 t +d 2 r --\^, (5.1) 

iTlJ-oo Jo \ r z ) 



where ip(r, t) is a solution to ( |2.15| ). Since the integrand in (|TTT|) is a total divergence we can 



also write Q(R) as a surface integral. Using fl2.15a| ) we have that 

oo R 



-oo o 

or 

,2 



«< fl > = ^/4*Mr^h (5 ' 2) 



9 f j i Sdrlp 



Q(R) = ( dt 

7T J 



-oo 



(5.3) 

r=R 



-2gS 

since the contributions from the t = ±oo surfaces and the r = surface vanish. 

We are interested in R much larger than any length scale associated with the solution. 
At this large value of R, the integral in (|5.3|) is only nonzero when the incoming pulse passes 
R at t ~ — R and when the outgoing pulse passes R at t ~ R. For t ~ —R we have 
ip{R,t) = ip p (R + t) + 0(1/ R) and 6{R,t) = 5 P (R + t) + 0(1/ R) while for t ~ R we have 
ip(R, t) =ip f (t-R) + Oil/ R) and 6{R, t) = 8 f (t - R) + 0(1/ R). Here ip p (u) = ip{u, -oo) 
and S p (u) = S(u, —oo) are the incoming pulse functions described in section |TTT| and ij)f(v) = 
ip(oo,v) and Sf(v) = 5(oc,v) are the associated future pulse functions. Substituting into 
and taking the limit as R goes to infinity gives 



Q = l 7 dv SmMM) _ l J d JWMM } . (5.4) 



7T 



1 — 2g5j(v) J Ti J [ 1 — 2g5 p (u 

Note that topological charge develops only to the extent that the outgoing profiles differ 
from the incoming profiles. 

We can develop ( |5.4| ) perturbatively in g if we use the perturbative expansion for the 
solutions ip(u,v) and 6(u,v). To zeroth order in g, ip(u,v) and 6(u, v ) are given by (|4.3a|) 
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and ( |4.3b|) so we see that to this order, the future and past profile functions agree, that is, 
il)fo(x) = ip p o(x) and Sf (x) = S p o(x). This means that the topological charge vanishes to 
second order in g. 

If we look at ( |1.2| ) along with ( |1.3| ) and ( |1.4j ) we see that the second order topological 
charge, outside the spherical ansatz, can be written as 

Q l2] = if^a E / d ' x ^ aP d^d a A a p , (5.5) 

where 

d»(d,A«-d u A;)=0 . (5.6) 

Thus the second order non-Abelian topological charge is equivalent to the topological charge 
in an Abelian theory. In appendix A we show, without restricting to the spherical ansatz, 
that for localized energy pulses which solve Maxwell's equations, the topological charge 
vanishes. 

To get the third order topological charge we need ipi(u,v) and 5\(u,v) whose general 
form is given by ( f4.10[ ) plus possible homogeneous pieces. With ipi and 8% given by ( (4.10[ ), 
we see that both go to zero for large positive u and for large negative v so that ipfi — ip p \ = 
Sfi = 5 p i = 0. Any homogeneous piece in ip\ and 5\ will make cancelling contributions to 
( |5.4|) so we have that to order g 3 the topological charge vanishes within the spherical ansatz. 

The fourth order topological charge is more difficult to evaluate. For arbitrary initial 
profiles we would need to find ip2 and $2 which entails evaluating ( |4.8|) with the sources ( |4.6c|) 
and ( |4.6d|) where ipo an d #o are given by ( [4.3|) and ipi and 5\ are given by ( |4.10| ) plus possible 
homogeneous terms. We have done this only in examples. For the examples considered we 
find that ^2 and 5/2 are non- vanishing even though ifj p2 = 5 P 2 = 0. Nonetheless, in every 
case the fourth order topological charge vanishes because of cancellations in the first term in 
( p.4| ). Our goal is not to prove that the fourth order topological charge is zero for arbitrary 
initial profiles. Our goal is to show that there are initial profiles for which Q is of order g n 
for some finite n. Therefore, onward to the g 5 term. 
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We have evaluated the fifth order topological charge for several initial profiles. In every 
case, we have found a nonzero result. In principle, once the initial conditions are specified 
the calculation of the fifth order topological charge follows from (|5.4|) and the results of 
section IV. In practice, at this order the calculations are rather involved and we relegate the 
details to an appendix. Here, we will describe the initial conditions we use, and quote the 
resulting for two examples. 

Initial conditions are specified by the functions 7 P (w) and a p (u) of equations Q3.9|) which 
we require to be g- independent. We make the choice 7 P (w) = because it simplifies the 
calculation considerably. From ( j3.9|) and ( }4.3| ), this choice implies G'{ ) and 

H(x) = 0. Because H(x) = 0, equation ( j4.10a| ) gives ipi{u,v) = 0, which agrees with the 
order g term of ( |3.9b| ). In general, the equations of motion ( p,15|) are invariant under the 
transformation g — > —g and 5 — > —5 with ip fixed. For 7 p = 0, the initial conditions ( |3.9| ) 
are invariant also. This means that for 7 p = 0, the expansion in g of 6 can contain only odd 
powers of g, while that of ip can contain only even powers of g. The order g term of ( |3.9a|) 



requires that S p i(u) 



1 G'(u) 2 . Therefore, we must add the homogeneous term 



6^ om (u,v) 
to 5 ± of ( ||T0g ) to obtain 

5i(u,v) = 



G\uf + G'(vf 



u — V 



dxG\xf 



(5.7) 



llG'iu)-^)?^-^ 2 



[u — V) 



■[G'{u) + G'{v)}[G{u)-G{v)} + 



dx G'(x) 



(5.8) 



it — V ' u — V 

Thus, choosing 7 P = means that S and ipx are zero f° r & H u an d v i while ijj (u,v) and 
5i(u,v) are given in terms of the single function G(x) by equations ( 4.3a ) and ( |5.8| ). 

In appendix B, we show that for 7 p = 0, the fifth order term in the topological charge 
( |5~4D can be written as the double integral 

6 



Q[5] = £! f°° dudvl- 
^ 71 J-oo [2 



-d^ d?8 X + 

2 [u — v) z 



2G'(u)G'(v) + 
4 



:[G(u)-G(v)Y 



3 ipoSx 



(u — v) 2 
1 V^o 5 1 



(u — v) 2 
[G'(u) + G'(v)][G(u)-G(v)} 



2 (u — v ) 2 (u — v)'" 



(u — v)' 



(5.9) 



We have evaluated equation ( |5.9|) analytically for several functions G(x). Here, we present 
the results for two examples. For 

G(x) = - - X 9 . (5.10) 



we find 



For 



l + x 2 



= -—£. (5.H) 
^ 1024 7T v ; 



; x < -1 

G"(s) = <( (1 - x 2 ) 2 ; -1<x<1. (5.12) 
; 1 < x 



we find 



[5] = 1336082432 £ (5 ^ 

^ 11402015625 vr v ' ; 

We have demonstrated that spherical solutions to the Minkowski space classical equations 
of motion for pure SU (2) which propagate freely in the far past can have nonzero topological 
charge. The power series expansion of Q in g can be nonzero at finite order — in the examples 
we have presented, at order g 5 . 

VI. FERMION PRODUCTION IN A BACKGROUND FIELD WITH 
FRACTIONAL TOPOLOGICAL CHARGE 

The solutions which we consider in this paper have the property that the topological 
charge inside a sphere of radius R, Q(R) given by (|1.1|), approaches a fixed value as R gets 
large and the limiting value, Q = lim^ 00 Q(i?), need not be an integer. In fact we have 
shown that Q can be developed as a power series in g so that Q can be made small by 
making g small. Suppose that we have a quantized left-handed fermi field, coupled to 
this classical non-Abelian background. Without introducing a Higgs field, we cannot give 
the fermion a gauge invariant mass, so we take the fermion to be massless. Due to the 
anomaly, the fermion number current J M = : vl^^ : is not conserved, that is, 
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d»> = J^e^Tr(F^) . (6.1) 

The question arises of how many fermions are produced in the presence of a background 
field of the kind we have described. 

Define the fermion number operator, inside the sphere of radius R, as 

N(R, t) = J Q R r 2 dr J d cos 9 d(j> J°(r, 9, (j>, t) (6.2a) 

and the flux operator, through the two-sphere of radius R, as 

F(R, t) = JdcosO d(j) R 2 J r (R, 6, 0, t) . (6.2b) 

(We will for the moment ignore the fact that these are not well defined operators and return 
to this point later.) Now, integrating (|6.1| ) from t = — T to t = T and over the inside of 
the sphere of radius R, we obtain 

N(R, T ) - N(R, -T ) + F dt F(R, t) = -f^ dt [ d 3 x Tr {F, u F a(3 ) . (6.3) 

J-T 6Zn- J-T J\x\<R 

We are interested in background shells which at early times are moving along r = —t, enter 
the sphere of radius R at t ~ — R, leave at t ~ R, and then continue to move outward 



along r = t. If we take T much larger than i?, then the right hand side of (|6.3| ) is well 
approximated by Q(R) of ( |5.1D and we can write 



N(R,T ) - N(R, -T ) + [ T ° dtF(R,t) = Q(R) . (6.4) 

J-T 



We are interested in tracing the evolution of a normalized fermion state |s) in the back- 
ground of a shell which enters and leaves the sphere of radius R. Pick the state \s) such 
that at t = —T there are no fermions inside the sphere, which implies 

(s\N(R,-T )\s) = 0. (6.5) 

The non-Abelian shell may have fermions riding along with it which enter and leave the 
sphere of radius R making (s\F(R, t)\s) nonzero for t ~ —R and t ~ R. Any fermions 
anomalously produced within the sphere will also contribute to (s\F(R,t)\s) when they 
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leave the sphere. For t ^> R, the energy shell is far outside the sphere and the background 
field is once again pure gauge within the sphere. Therefore, at late times fermions inside 
the sphere are free, and if we wait long enough they will leave the sphere. Thus, for T 
sufficiently large, we have 

(s\N(R,T )\s) = 0. (6.6) 

Taking the expectation of (|6.4j) in the state \s), we obtain 

[ T ° dt(s\F(R,t)\s) = Q(R) • (6.7) 

We interpret ( |6.7|) as follows. Imagine a spherical detector corresponding to the operator 
( |6.2b|) , which at any time t measures the flux of fermions passing out through the two sphere 
of radius R. (Fermions going in make a negative contribution, as do anti-fermions going out.) 
The detector is in continuous operation from t = —Tq to t = Tq and therefore measures the 
total fermion number passing through the two sphere. For a given background field, we 
cannot predict the outcome of an individual measurement of the total fermion number 
passing through the two sphere. However, if we repeat the measurement many times with 
the same background field and the same state \s), then the average of the observations will 
be Q(R). 

Although we believe that the spirit of the argument just given is correct, it is specious 
because the operators N(R,t) and F(R,t) have infinite vacuum fluctuations, that is, 

(0|iV 2 (it:,f)|0) = oo (6.8a) 

and 

(0\F 2 (R 7 t)\0) = oo . (6.8b) 

Therefore, these operators also have infinite fluctuations in any state in the Fock space; in 
particular (s\N 2 (R, t)\s) is infinite. Thus, we cannot sensibly discuss the measurement of 
these operators. 
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The infinite vacuum fluctuations arise because of the sharp cutoffs in the r integration 
at r = R and in the t integration at t — ±T in expressions ( |6.2a| ), ( |6.2b| ), and (|6.3|). This 
can be remedied by introducing a smoothing function. Let 

J y ' \ 0, if r > R D or \t\ > T D v ' 

and let f{r,t) fall smoothly from 1 to as r varies from R to Rd and as \t\ varies from To 
to Tjj. Define the smoothed operators 

N f = - r dt r r 2 dr f d cos 9 d<f) ^jf' ^ J°(r, 0, 0, t) (6.10a) 

J —oc JO J C/t 

and 

T) = dtj\ 2 dr J d cos 6 d<f) Q J r (r, 0, 0, t) . (6.10b) 

Because df(r,t)/dt is only nonzero for T < |i| < To, it is also convenient to define 

A> f early = /~ T °rft r r 2 dr f d cos dxj> ^ ^ J°(r, 0, 0, t) (6.11a) 



7 -J_ Td ^ Jq > u,j 
and 

A>' ate = - /^dt fVdr / dcos6d<P^^j°(r,6,<f),t) (6.11b) 
^ JT jo j at 

where iV^* 6 - iV / eaxly = iV>. Note that the eigenvalues of N f , F f , iV / early , and iV^ 6 , are not 
integers. A short calculation shows that the vacuum fluctuations of the smoothed currents 
are given by 

d 4 x f d*yf(x)f(y)(0\>(x)r(y)\0) = ~ / 0(k°) 9(k 2 ) \f(k)\ 2 (Ff - gftf 



(6.12) 



where / is the fourier transform of /, and the currents are understood to be normal ordered. 
From ( |6.12| ), we see that by a judicious choice of f(r,t), the vacuum fluctuations can be 
made arbitrarily small. For the massless fermion case which we are considering this requires 
Rd to be much larger than R and To to be much larger than To. We also require that 
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T ^> Ro so that during the measurements of A^ carly and Nj atc the energy shell is far 
outside the detector. To summarize the relations between the relevant length scales, we 
have A <C R <C R^, <C T <C Tp where again A is the length which characterizes the pulse. 
If we multiply both sides of (|6.1|) by f(r,t) and integrate, we get 

^late _ pearly + p f = jt_ ^ £° ^ J fa f ^ f) Tr ( F ^ F ^) . ( 6 . 13 ) 

Using the equations of motion for the fields in the spherical ansatz, we can write the right 
hand side of fl6.13| ) as 

f r dt f°° dr f(r, t) (-d 2 t +dl--\^. (6.14) 

/7T J —oo JO \ T J 



For r ^> A, and in particular for r > R, the integrand of Q6.14 ) is well approximated by 



f(r, t) {—d^ + df) i/j which is near zero because wherever and whenever ip is nonzero at large 
r, it is given either by ijj p (t + r) or by ip/(t — r) up to terms of order 1/R. Also, for \t\ > T 
the shell is outside r = R and so at these early and late times ip is close to zero for r < R. 
Thus, the integrand of (|6.14 ) is close to zero wherever and whenever f(r,t) ^ 1. Therefore, 



the expression ( |6.14|) is Q up to terms of order 1/R. 



Taking the expectation of ( |6.13|) in a state \s) which has 

(s|A>/ arly |s) = (6.15) 

and neglecting terms of order 1/R, we obtain 

(s\F f \s) = Q, (6.16) 

where again we argue that (s|A^j ate |s) = for sufficiently large To. As before, we interpret 
( |6.16 ) in an expectation value sense. However in this case if the fluctuations of the flux 



operator are small enough, then each observation of the net number of fermions flowing 
through the detector corresponding to the smoothed operator ( 6.10b ) will yield a number 



close to Q. In this sense, fermion number violation is seen with each observation. 
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VII. CONCLUDING REMARKS AND OPEN QUESTIONS 



We have shown that spherical solutions to the classical equations of motion of SU(2) 
gauge theory which propagate freely in the far past have nonzero topological charge. Fur- 
thermore, the topological charge of such solutions can be developed as a power series ex- 
pansion in g. We presented examples for which Q has a nonzero contribution at order g 5 . 
In these examples Q is nonzero for ^-independent initial conditions. In particular, this sug- 
gests that in a corresponding quantum process, the number of particles in the initial state 
is ^-independent, and need not be large at small g. We showed that if we couple quantized 
chiral fermions to classical gauge field solutions with nonzero Q, anomalous fermion number 
violation will occur. In a series of measurements using the operators of section VI with the 
same background gauge field and the same quantum mechanical fermion state, the average 
number of fermions anomalously produced is given by the topological charge Q. For the 
solutions we have discussed, the topological charge vanishes as a power of g. This suggests 
that anomalous fermion number violation in high energy scattering processes may occur at 
finite order in g, rather than being exponentially suppressed. 

We cannot yet apply our results directly to scattering of W bosons at energies much 
higher than the sphaleron energy because a number of important open questions remain to 
be considered. These questions arise because we have done a classical calculation in the 
spherical ansatz for an SU{2) gauge theory without a Higgs field. Let us consider each of 
these idealizations in turn. 

First, in everything we have done, we have treated the gauge field classically, whereas we 
are ultimately interested in quantum scattering processes with a two-particle initial state. 
We need to relate our classical results to quantum mechanical scattering. One possible 
approach is to relate small amplitude classical solutions to quantum mechanical coherent 
states with small mean particle number. In this context, it is interesting to note that the 
^-expansion we have done is equivalent to an expansion in the amplitude of the initial 
pulse. Therefore, the topological charge at order g 5 is proportional to the fifth power of the 
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amplitude of the initial pulse, and is not exponentially small for small amplitude pulses. 
This suggests that there may be scattering processes with few particles in the initial state in 
which Q is not exponentially suppressed. Another way of addressing the problem might be 
to relate the lowest order term in a semi-classical expansion of the ^-matrix for a quantum 
scattering process to the solutions of classical equations of motion for a specified initial 
configuration. The quantum lessons which our classical results are trying to teach us remain 
to be learned. 

Second, we have worked within the spherical ansatz. This is a reasonable starting point 
since it makes the calculation tractable, and typically the lowest partial waves dominate 
scattering processes. An important open question here is that of stability. We have seen 
no evidence that the solutions we discuss are unstable to perturbations within the spherical 
ansatz. However, we have not checked their stability to non-spherical perturbations. 

Finally, we have worked in SU (2) gauge theory without a Higgs field. It remains to be 
determined how the inclusion of the Higgs field affects our results. It is possible that in 
scattering at very high energies the W mass and sphaleron barrier are irrelevant and our 
results are a good approximation to results in the full theory. However, including the Higgs 
field changes the theory in qualitative ways. The theory is no longer scale invariant. The 
behavior we describe at early and late times — pulses propagating without distortion - 
will not occur once the Higgs field is included. Thus, while our results are likely applicable 
in the full theory in the limit of infinite energy scattering, it remains to be seen how they 
can be applied at high but finite energy. 
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APPENDIX A: ABELIAN TOPOLOGICAL CHARGE 



Here we show that for a localized solution to Maxwell's equations in empty space, the 
topological charge vanishes. We do not use the spherical ansatz. It is convenient to work 
with the gauge invariant electric and magnetic fields E(x, t) and B(x, t) which obey 



V-E = 
V-B = 

The most general solution can be expressed as 



* ft dE 



(Ala) 
(Alb) 



E(x, t) = Jd*k [eik) exp (iut - ik • x ) + e *(k) exp + ik 



x 



(A2a) 



and 



B(x,t) 



d 3 k r 



k x e(k) exp (iuut — ik ■ xj + k x e *(&) exp (^—iut + z'A; ■ 



(A2b) 



where u; = \k\ and e*(fc) is an arbitrary function of k subject to k ■ e(k) = 0. 
Now the Abelian topological charge can be defined as 

1 



Q 



Abelian 



(fx E-B 



(A3) 



which in terms of E and B given above is seen to be 

oo 

gAbclian = f „ fa f d COS 6 difd(2uj) [?(£) • (jfe X ?(-£)) + C.C. 



(A4) 



The delta function in u comes from the / dt in ( |A3| ) . Because of this delta function, the 
integral will vanish unless e(k) goes like 1/w as w -> 0. If we know E and B at one time, 
say t = 0, we can determine e*(fc): 



*r(*) 



(i 3 x 
2 J (2^)3 



£(£ , 0) x B(x, 0) 



exp ik ■ x . 



(A5) 



As goes to zero, e(k) will not blow up if 



d 3 x E(x, 0) < oo 



(A6a) 
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and 

J d 3 x B(x, 0) < oo . (A6b) 

For localized pulses these conditions are satisfied and the Abelian topological charge van- 
ishes. Note however that finite energy does not ensure ( |X6| ) so vanishing Abelian topological 
charge is not guaranteed under the most general circumstances. 

APPENDIX B: EVALUATING THE FIFTH ORDER TOPOLOGICAL CHARGE 



In this appendix we derive equation ( |5.9|) for the fifth order topological charge for initial 
profiles for which 7 P = and a p is independent of g. As was explained in the text, this implies 
that ip n = for n odd and 5 n = for n even. We first need the third order contribution to 
the integrands of equation 



From (|5.8|) we see that as u — > oo, 5\ — > — (l/2)G"(t>) 2 and as v — > — oo, 8\ — > —{l/2)G'{u) 2 . 
Therefore, the last term in makes cancelling contributions to the two terms of Q5.4D . 
Furthermore, there is no contribution from the first two terms of ( [Bl|) asw^ — oo since 63 
and ip2 vanish in this limit. Thus, 



q[5] = t_ r dv (-^ d v s 3 + s*°™dM 

71 J-00 



m=oo 



(B2) 



where the first term has been integrated by parts, and in the second term we have used the 
fact that at large u we can replace S\ by the homogeneous solution 5i° m given by ( |5.7|) . 
These single integrals may be turned into double integrals using 

/OO /"OO rU 

dvr)d v & — du dvrjl , (B3) 
-00 u=oo J— OO J— OO 

where $, given by ( |4.8| ), is the solution to equation ( |4.7| ) with source J, and rj is any 
homogeneous solution of ( [4.7] ). This relation can be established by substituting ( |4.8| ) into 
the left hand side of (|B3| ) and writing the homogeneous solution 77 as in equation (|4.3| ). 
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Using (p3|) we write (B2) as 

)[5] _ 9_ 



/oo pu 
du dv(- 4j 4 + Si° m lf) 
-oo J — oo 



7T 



(B4) 



Using the equations of motion, ( |4.6c| ) can be written 



(u — t>) ; 



(B5) 



Similarly, 



J 3 5 = -cU - -(Ji d^ - -d^ d^ 2 + 



46 2 



(B6) 



(u — v) 2 

Using equations ( |4.3a| ), (|5.7|), ( |5.8| ), (P^|), and ( |BT)| ) we can evaluate every piece of (B4) 
except the part which has the term from flB6|) involving ip 2 . Integrating this term by parts 
yields 



1 



oo J — oo 



oo ru 



oo J — oo 



du / dv d^ipo d^ijj 2 = du dv ip 2 d u d v % - - \ dv tp 2 d v ip Q 



(B7) 



Using ( |5.7| ) and (|B"3|), the surface term can be written 



\r 

£ J —oo 



dv i\) 2 d v tpQ 



dv ip 2 d v 5l 



horn 



/oo ru 
du / dv 5\ om lt . (B8) 
-oo J —oo 



Note that this term cancels the second term of (|B4|). We now examine the first term on the 
right hand side of (p7|). Using ( |4.8|) we find 

ru ru rv 

du I dv I du' I dv' F(u,v;u',v')I 2 (u',v') 

d u d v ip 2 (u,v) 

du' I dv' I du I dv F(u,v;u',v')d u d v ipQ(u,v) 



oo J — oo 



du I dv ip 2 d u d v ip 2 



—oo J —oo Jv 



oo J —oo 



oo ru 



ii(u',v'). 



(B9) 



We perform the inner u-v integrals as follows. Define the function 

/oo ru' 
du I dv F(u,v;u' ,v')d u d v ip 2 (u,v) 

One can show that 



(BIO) 
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d u d v + 



(u — vy 



(Bll) 



Hence Q(u,v) is a solution to the inhomogeneous equation ( [Blip with the boundary condi- 
tion 



lim Q(u, v ) = 



(B12) 



Note that the right hand side of ( [Blip is the same as the right hand side of (|4.6a|) with 
= —i>o- Also, the solution ipi given by (|4.10a|) satisfies the condition ( B12|) . Therefore, 
<d(u,v) is given by (|4.10a|) with H = —G: 

6 



G(u,v) = 2G\u)G\v) + 



(u — V) 



■[G(u)-G(v)Y 



(u — v) : 



: [G'(u) + G'(v)][G(u)-G(v)] . 



(B13) 



Thus, 



du 

oo J— CO 



u poo ru' 

dv ip2d u d v ipl = du' dv' e(u',v')lf(u',v') . 

-OO J —CO J—oo 



(B14) 



Using (|B~6|), (P?D, (|B8|) , ( |B13| ), and ( |B14| ) in ( P4| ) and integrating by parts several times we 
find 



Q 



dudv 



7T J-oo 



1 n r V 3 t^l 



1 t 



(B15) 



.2~~ z 2(m-w) 2 (u-v) 2 8 (u-v) 2 . 

which is equation ( |5.9|) . In (P15|) we found it convenient to use the fact that the integrand is 
symmetric under the interchange of u and v to extend the integration range over the entire 
u-v plane. With equation ( P15|) in hand, we can choose G(x), obtain ip Q from ( |4.3a|) , Si 
from ( |5.8|) , if from (|B5|), and from ( P13|) , and then evaluate the fifth order topological 
charge. 
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